Charge degrees in the quarter-filled checkerboard lattice 
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For a systematic study of charge degrees of freedom in lattices with geometric frustration, we consider spin- 
less fermions on the checkerboard lattice with nearest-neighbor hopping t and nearest-neighbor repulsion V at 
quarter-filling. An effective Hamiltonian for the limit \t\ <g V is given to lowest non- vanishing order by the 
ring exchange (~ t^/V^). We show that the system can equivalently be described by hard-core bosons and map 
the model to a confining (7(1) lattice gauge theory. 

PACS numbers: 



Recently, there has been enormous interest in the proper- 
ties of quantum magnets on lattices with frustrated geometries 
where the number of classical ground state configurations in- 
creases exponentially with the number of sites. The arguably 
most spectacular phenomenon in this field is fractionalization 
of quantum numbers (e.g., spinons), see Refs. IllQl and cita- 
tions therein. Now, attention shifts to charge degrees of free- 
dom. A model on the checkerboard lattice has been intro- 
duced in which an added fermion can decay into two mobile 
pieces which carry the fractional charge of e/2 each jst]. Nu- 
merical studies of this model have shown that the fractional 
charges are weakly confined and form bound pairs of large 
spatial extent. These are expected to lead to interesting physi- 
cal effects . In the present work, a related model, namely 
the checkerboard lattice on which one-quarter of the sites are 
occupied by strongly interacting spinless fermions, is consid- 
ered. We dicuss an effective Hamiltonian that describes the 
low-energy physics. This Hamiltonian has a U{1) gauge in- 
variance which we use to relate the problem at hand to the 
compact quantum electrodynamics in 2 + 1 dimensions (0]. 

Our starting point is the Hamiltonian of spinless fermions 
in the checkerboard lattice with nearest-neighbor hopping am- 
plitude t and a strong repulsive interaction V on neighboring 
sites. The low-energy manifold at quarter-filling in the limit of 
strong interactions (|t| ^ V^) is given by those configurations 
which have exactly one fermions on each crisscrossed square. 
This local constaint corresponds to the so-called tetrahedron 
rule, which applies in the half-filled case IH ItIi. The low- 
energy manifold of states fullfilling the local constraint can 
alternatively be represented by a model on the square lattice 
connecting the centers of the crisscrossed squares. The parti- 
cles sit in the middle of the dimers (see Fig.[T|(a)). The condi- 
tion of one fermions per crisscrossed square is being translated 
into the requirement of one dimer touching each site. Hence, 
states satisfying this constraint are represented by hard-core 
dimers on the square lattice. Notice that this model is very 
reminiscent of quantum dimer models on the square lattice 

i. 

In lowest non-vanishing order in t/V, the Hamiltonian for 
spinless fermions on a quarter-filled checkerboard lattice can 
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FIG. 1: (a) Square lattice connecting the centers of the crisscrossed 
squares of the checkerboard lattice with links. One exchange pro- 
cess around a hexagon is shown explicitely. (b) Labeling of links on 
flipable double-plaquettes at position x with unit vectors ei and e2. 



be reduced to the following effective Hamiltonian 10]: 
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The sums are performed over all polygons of perimeter six, 
and g ^ /V'^. The dots indicate particle locations. There 
is a difference with respect to the 'bosonic' quantum dimer 
model lUl] resulting from the fermionic statistics: The small- 
est resonating plaquette has perimeter six, rather than four It 
is noticeable that no Fermi sign problem remains. Only those 
processes that have the same overall sign will occur in the 
ground state of the quarter-filled case In passing we men- 
tion that the overall sign of the ring-exchange amplitude g in 
([T]i can be chosen arbitrarily. ifTlll 

The effective Hamiltonian ([T|) conserves the number of 
fermions on each crisscrossed square. Consequently, the num- 
ber of dimers touching each site of the lattice connecting the 
centers of the crisscrossed squares is not changed. This con- 
servation generates a [/(I) gauge invariance, as it is usually 
the case in dimer models \%. The gauge structure suggests 
that we can gain further insight with respect to low-energy 
excitations by writing our model as an U{1) lattice gauge the- 
ory. The usefulness of this approach has already been shown 
for the quantum dimer model (QDM) on the square lattice |0] 
as well as for three dimensional spin systems jlOll . 

Define for each link (x, x + Cj ) of the square lattice an in- 
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teger variable rij (x), where x denotes the coordinates of a lat- 
tice site and ej=i 2 are unit vectors along the axes as shown in 
Fig. [lib). The states |{nj(x)}) span an enlarged Hilbert space 
which has integer numbers for the links instead only zero or 
one. We can consider the states \{nj{x)}) as eigenstates of 
quantum rotor operators fij (x) with eigenvalues Uj (x). 

In order to express the effective Hamiltonian in terms of 
the rij(x), we introduce phases 0j(x) e [0, tt) on the links 
which are canonical conjugate to 7ij(x). Using the fact that 
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act as ladder operators, we can write 
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Here, the argument of the cosine term contains the sum over 
phases with alternating signs around the polygons of 

perimeter six (double plaquettes). In the limit U/g 00 
Eq. (|2]i is a faithful representation of the effective Hamiltonian 
O- 

We introduce staggered gauge and electric fields on the bi- 
partite square lattice x = (xi, X2) by 



shown in Fig.[T] Eq. (|4| has similarities with the Hamiltonian 
of the compact quantum electrodynamic (QED) in 2 + 1 di- 
mensions in which the considered charges correspond to frac- 
tional charges of e/2 i^. Polyakov showed for the compact 
QED in 2 + 1 dimension that it has an unique and gapped 
ground state. Two charges are confined and the energy grows 
linearly with the distance between the two charges. Our model 
shows important differences. The fields Ej{x.) are half inte- 
gers instead of integers and the constraint selects configura- 
tions with a background charge p(x). This leads to a frustra- 
tion which is reflected by the macroscopic degeneracy of the 
classical ground states. Furthermore, the definition the flux in 
the cosine term differs. 



Formulation ^ provides an excellent starting point for fur- 
ther systematic investigations. The plaquette duality transfor- 
mation allows to map the Hamiltonian to a height model and 
to use path integrals for a detailed study of the ground-state as 
well as low-energy excitations |@]. 



^,(x) = (-l)-^+-^0,(x) 
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The local constraint that each site is touched by exactly one 
dimer reads 

(a,4 (x) - p(x)) iPhys.) = 0, p(x) = (-l)(-^+--^+i),(3) 
where the lattice divergence is defined as 

A,E, (x) = ^i(x) -Ei{x- ei) + £2(x) - E^i^ - 63). 

The constraint is now reflected by the standard Gauss' law ^ 
in the presence of a staggered background charge density. The 
Hamiltonian (|2]i in staggered variable reads 



n^ff=Uj2E]ix)-2gJ2cos 
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Here, the argument of the cosine term denotes the oriented 
sum of staggered vector potentials Aj (x) around double pla- 
quettes. The vector x labels the site in the left bottom as 
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